The concurrent work G2Gs presents a similar two-step framework, while our method is more general and scalable. We keep muted about G2Gs before its conference version is available since we have some concerns about it. This discussion will be included into our revised version. Atom mapping is optional for our method. The synthon approach can also work for reactions without provided atom mapping (L208-212).

We thank all the reviewers for insightful comments and suggestions. We address the two remarks below. This enables the adaptive minimaxity for Sobolev and Holder classes. Our results do not directly apply to that stronger setting. Thanks for the detailed and insightful review.